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Abstract 

We consider a class of non-linear PDE systems, whose equations pos- 
sess Noether identities (the equations are redundant), including non- varia- 
tional systems (not coming from Lagrangian field theories), where Noether 
identities and infinitesimal gauge transformations need not be in bijection. 
We also include theories with higher stage Noether identities, known as 
higher gauge theories (if they are variational). Some of these systems are 
known to exhibit linearization instabilities: there exist exact background 
solutions about which a linearized solution is extendable to a family of ex- 
act solutions only if some non-linear obstruction functionals vanish. We 
give a general, geometric classification of a class of these linearization 
obstructions, which includes as special cases all known ones for relativis- 
tic field theories (vacuum Einstein, Yang-Mills, classical N — 1 super- 
gravity, etc.). Our classification shows that obstructions arise due to the 
simultaneous presence of rigid cosymmetries (generalized Killing condi- 
tion) and non-trivial de Rham cohomology classes (spacetime topology). 
The classification relies on a careful analysis of the cohomologies of the 
on-shell Noether complex (consistent deformations), adjoint Noether com- 
plex (rigid cosymmetries) and variational bicomplex (conserved currents). 
An intermediate result also gives a criterion for identifying non-linearities 
that do not lead to linearization instabilities. 

1 Introduction 

It is well known that, on spatially compact manifolds, the solution space of 
Lagrangian gauge theories like General Relativity (GR) and Yang-Mills (YM) 
theory is an infinite dimensional manifold with singularities [T| [5] . The same 
phenomenon can occur in the solution space of other non-linear systems of 
partial differential equations (PDE systems), be they Lagrangian field theories 
or not. A background solution is said to be linearization stable if a solution space 
neighborhood of it can be modeled on the vector space of linearized solutions 
about that background. That is, for every linearized solution, there exists a 
1-parameter family of exact solutions tangent to it. Otherwise, the background 
is said to be linearization unstable, in which case there exist linearized solutions 
not tangent to any 1-parameter family of exact solutions. 



In the case of GR, the solutions space singularities at linearization unstable 
backgrounds are of conical type. That is, the corresponding solution space 
neighborhood can be modeled on the zero set of a quadratic constraint on the 
space of linearized solutions. The presence of such constraints is intimately 
linked with the spatial compactness of the underlying manifold. For spatially 
non-compact manifolds with asymptotically flat boundary conditions, conical 
singularities are absent [3J 2] . 

Based on a remarkable observation of Moncrief [SJE], the necessary and suffi- 
cient conditions for for linearization stability (resp. instability) have been identi- 
fied with the absence (resp. presence) of Killing vectors on the background [2 [2]. 
Similar results have been obtained for other gauge theories, including Yang- 
Mills [7], Einstein- Yang-Mills [8] and supergravity [S]. In each case, the notion 
of a Killing vector needed to be generalized to express the sufficient conditions 
for linearization instability. 

In this note, we generalize the Killing condition to a general class of La- 
grangian gauge theories, and even to non-Lagrangian theories, that satisfy non- 
trivial Noether identities [101 111) (a.k.a gauge identities |12) or compatibility 
identities |13l 114] . which mean that the field equations are redundant). In La- 
grangian gauge theories, the Noether identities precisely dual to generators of 
gauge transformations, by Noether's theorem. In non-Lagrangian theories, the 
two need not be connected [11] . We then show how these theories acquire poten- 
tial linearization instabilities, as in the above examples. In fact, this shows that 
it is the Noether identities that are responsible, rather than the gauge symme- 
tries. Thus, since the classical Killing condition is associated with symmetries, 
we call the generalized condition responsible for linearization instabilities co- 
Killing. Since, as already mentioned, Noether's second theorem links gauge 
symmetries and Noether identities, this distinction may have been ambiguous 
in the past. Briefly, fields satisfying the co-Killing condition are those in the 
kernel of a adjoint linearized Noether operator and are called co symmetries [15] . 
In specific cases, they have also been called reducibility parameters [12] . 

Starting with the original work of Fischer & Marsden on vacuum GR |16| . 
the conditions for linearization instability have often been detected through an 
analysis of the constraints imposed by the field equations on the initial data. 
In 4-dimensional GR, this requires the well-known 3+1 ADM decomposition. 
However, given that the resulting conditions (presence or absence of Killing 
vectors) are spacetime covariant it is desirable to have a covariant derivation of 
the standard and generalized co-Killing conditions as well. Such derivations are 
indeed available for specific examples as for instance for vacuum GR in [1] (see 
also [T7], where similar ideas appear outside the specialized literature). In this 
paper, we give a covariant derivation of the generalized co-Killing condition, 
applicable to the same large class of theories mentioned above. 

Since we analyze a large class of field theories instead of working on specific 
examples, we require an adequate abstract framework to state the necessary 
hypotheses and carry out the analysis. The abstract framework is provided by 
the jet-bundle formalism and various associated differential complexes and their 
cohomologies. One benefit of the abstraction is to attract attention to the fact 
that the necessary hypotheses are surprisingly liberal. In particular, the class 
of admissible PDE systems is larger than the determined elliptic and hyperbolic 
systems or gauge theories closely resembling those. Also, it becomes clear how 
the topology of the spacetime manifold gives rise to linearization instabilities 
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and that any non-trivial de Rham cohomology class may be responsible (not 
only one dual to a compact Cauchy surface). 

The jet-bundle formalism is introduced in Section^ its contents are standard 
and serve mostly to fix notation. The main hypotheses imposed on the PDE 
systems under consideration are described in Section 12.21 Section [3] introduces 
three important concepts: consistent deformations (Section I3.1J) . cosymmetries 
and the co-Killing condition (Section [3.21) . and null sources (Section [3. 3p . Sec- 
tion 2] contains the main results of the paper. Linearization instabilities and ob- 
structions are defined in 14. II Special conserved currents, deformation currents, 
are defined by pairing consistent deformations and null sources in Section l4~2l 
Our main result is obtained in Section 14.21 the deformation current defined by 
the non-linearity gives rise to a linearization obstruction valued in the de Rham 
cohomology of the spacetime manifold. The relation of this result to previous 
work is briefly discussed in Section 14.41 while Section 14.51 discusses obstructions 
generated by non-trivial asymptotic boundary conditions rather than non-trivial 
topology. Finally, the general result is specialized to several examples, some of 
which are novel, in Section [5] and Section [5] concludes with a discussion. 

2 Local geometry of PDEs 

2.1 Jets, local forms and the variational bicomplex 

The natural setting for the local analysis of differential equations is that of jet 
bundles. The needed concepts and notation are briefly introduced below. More 
details can be found in the standard literature; see for instance [T5J Q33 HOI I2~T1 

Hii. 

Given a vector bundle F — > M over a connected n-dimensional smooth man- 
ifold M , the k-jet bundle J k F — > M is a vector bundle whose defining character- 
istic is that for any (possibly non-linear) differential operator /: T(F) — > T(F') 
of order k, there exists a canonical factorization f[u] = / o j k u for any section 
u: M — ;> F, where the k-jet prolongation j k : T(F) — > T(J k F) is composed with 
a smooth bundle map /: J k F — > F', which by a slight abuse of notation we de- 
note using the same symbol as the original differential operator. Composing the 
differential operator / with an Z-jet prolongation canonically defines a new differ- 
ential operator pif : J l+k F — > J l F' called its I -prolongation, j l f[u] — pif o j k u. 
Given a trivializable restriction Fjj —> U of F to a chart U C M with local 
coordinates (x l ) and fiber-adapted local coordinates (x l ,u a ), there is a corre- 
sponding adapted chart J k Fjj C J k F with adapted local coordinates {x l ,itf), 
where / = i\---ii runs through multi-indices of orders |/| = I = 0, . . . , k. 
In these coordinates, the fc-jet prolongation is given by j k u(x) — (x l ,diu a (x)), 
while the /-prolongation is given by pif[u](x) = (x l , di f b [u](x)), where f[u] (x) = 
(x l , f b [u](x)) in fiber-adapted local coordinates (x l ,v b ) on F'. For any / > k, 
discarding the information about all derivatives of order > k defines a natu- 
ral projection J l F — > J k F. The projective limit J°°F := jhn^^ J k F de- 
fines the oo-jet bundle. The oo-jet prolongation j°° and co-prolongation 
arc defined in the obvious way. By composing with the natural projection 
J°°F — > J k F, the differential operator / also canonically defines the smooth 
bundle map /: J°°F -> J k F 4 F' , which is again denoted by the same symbol 
/. Conversely, due to the projective limit construction, any smooth bundle map 
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/: J°°F — > F' can only depend on finitely many coordinates of its domain, 
which means that there exists a k > such that this bundle map canonically 
factors as /: J°°F — > J k F — > F' , with the smallest such k being the order of /. 

Denote by TM and T*M the tangent and cotangent bundles of M. Also, let 
A k M = f\ k T*M be the bundles of fc-forms. Denote by fl k = n k {M) = T(A k M) 
the spaces of differential forms, with ft* = fe il k . We call fi*(J°°F) the space 
of local variational forms on F. The de Rham differential on Q*(J°° F) canoni- 
cally splits into the sum d = dh + d v , where the respective horizontal and ver- 
tical differentials are individually nilpotent and anti-commutative, d^ = d 2 = 
and dhd v + d v dh = 0. The defining action of the horizontal differential on 
adapted local coordinates (x l ,Uj) is dhX 1 = dx l and d^uf = u%dx l ] then sim- 
ply d v = d - d h . Denote by Q h '°(F), Cl°' v (F) c fl*(J°°F) the subspaces of 
local horizontal and local vertical forms, which generate the entire space of lo- 
cal variational forms by wedge products. Hence, we have a natural bigrading 
n*(J°°F) = Cl h,v (F), where < h < n and < v < oo. There is a 

natural inclusion Q k (M) C Sl k,0 (F), via the pullback along the natural pro- 
jection 7TOO : J°°F — > M, where the image of the inclusion is said to consist of 
field-independent forms. 

The operators dh and d v together with the horizontal- vertical bigrading turns 
the space of local variational forms, augmented as shown below, into the varia- 
tional bicomplex of F [191 120) : 
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o ► n°> 2 -^-> n 1 ' 2 ■ ■ • it- 1 * 2 -^-> ir< 2 — U t 2 ► o 



We have abbreviated f2 fe = VL k {M) and fl h ' v = Cl h ' v (F). All arrows commute, 
which defines the Euler- Lagrange derivative 5el = I°d v , where / is the so-called 
interior Euler operator. The space T k = (J) a d v u a A is called the space 

of local functional k- forms is defined by the image of I and the variational 
differential S v is defined as the unique map commuting with the rest of the 
diagram. 

We can define the cohomologies, ff '"(dh) and H h ' v (d v ), of the horizontal 
and vertical differentials in the obvious way; they correspond to the cohomolo- 
gies of the corresponding parts of the above rows and columns. All the columns 
and rows are exact, with the exception of the underlined nodes, with the proviso 
that the bent 8 el arrow makes the following sequence exact also except at the 
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underlined nodes: 

o > n Q -o^...n n - 1 ^^n n '°^T 1 ^T 2 ^--- . (2) 

At the underlined nodes, the horizontal cohomology is completely characterized 
by the de Rham cohomology of M, H k >°(d h ) ^ H$ R (M). At the same nodes, 
the vertical cohomology consists of all field-independent forms, iJ fe,0 (d v ) = 

n k (M) c n k >°{F). 

Note that any local horizontal form to £ Q*'°(F) naturally defines a smooth 
bundle morphism ui : J°°F — > A*M as well as a differential operator to : T(F) — >• 
r(A*M), u[4>} = to o j°°tp for ip £ T(F), where we have slightly abused notation 
by denoting each of these naturally associated objects by the same symbol to. It 
is convenient to generalize this construction by replacing the bundle A*M — > M 
with an arbitrary vector bundle H — s- M. Any differential operator /: J°°F — » 
H shall also be termed a local section of H (also F -local if such precision is 
necessary). Denote the space of all such local sections by Y F (FF). Evidently 
the space of local horizontal forms is the same as the space of local sections of 
A*M, n*'°(F) ^ r F (A*Af). The pullback along the bundle projection F -> M 
induces a natural inclusion T(H) Tp(H) of field-independent local sections 
in the space of all local sections. Similarly, given two vector bundles G — > M 
and H — > M, we call a differential operator /: IV (G) — > T F (H) local (or 
F -local) when fty,^}} = f ° 0'°°^ X Poo£), where $ £ T(F), £ £ T F (G) 
and on the right-hand side /: J°°(F X H) —> H. If f[ip,Q is linear in its 
second argument, we write it as /[?/'; C]- Of course, any differential operator 
/: r(G) — ► T(H) naturally defines a field-independent local differential operator 
f:T F (G)^T F (H). 

Local sections and local differential operators are discussed in [T31 HI] as 
C-modules and C-differential operators. 

2.2 PDE submanifold 

In the jet bundle setting, a PDE system has a dual description [2~T1 fT5l [HI 123] . 
as a non-empty, smooth sub-bundle £°° C J°°(F) over M and as the zero- 
set of the oo-prolongation Pooe: J°°F — > J°°G of a smooth bundle morphism 
e: J°°F — > G, where the vector bundles F — s> M and G — > M are respectively 
referred to as the field bundle and the equation bundle. The requirement of 
having both descriptions is non-vacuous. It is possible that the zero set of Pa^e 
is not a submanifold (say if the Jacobian of e is not of constant rank) and it is 
possible that £°° is not the zero set of any smooth bundle map (see [22j Section 
7] for the necessary and sufficient topological condition on the normal bundle 
of the embedding £°° C J°°(F)). A section <j>: M — > F is a solution of the 
PDE system if its oo-jet prolongation is contained in the PDE submanifold, 
j°°<j>(M) C or equivalents if e[<j>] = 0. 

Clearly the smooth PDE sub-bundle J°°F D £°° — > M provides an intrinsic 
description of the PDE system, while the choice of the equation form, the bundle 
G — > M and the map e, are not unique. Suppose that the differential operator e 
is of order k and denote the natural projection ir k : J°°F — > J k F. Without loss 
of generality, we can presume that the zero set £ of e: J k F — ► G is a smooth 
sub-bundle of J k F and that it agrees with the projection of £°°, £ = ir k £°°. 
Such an equation form is said to be involutive. In the sequel we freely use an 
equation form e[ip] =0 or the submanifold £ C J k F to define a PDE system. 
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Define the space of on-shell local variational forms n*(£°°) to be the pull- 
back image of r2*(J°°_F) along the inclusion £°° C J°°F and d the de Rham 
differential on it. The horizontal- vertical bigrading and the split of the de Rham 
differential, d = dh + d v , commute with the pullback, which immediately defines 
the decomposition f2*(£°°) = J^h » Q h ' v (£) and the on-shell variational bicom- 
plex with the projected operators dh, d v , Sel and 6 V . The cohomology Hg'°(db) 
in the space of on-shell local horizontal forms is also called the characteristic 
cohomology of the PDE system and is denoted #* har (£) QH |2TJ 125 ■ We sti11 
have a natural inclusion of the de Rham cohomology of M in the characteristic 
cohomology, if| R (M) C H* hal (£), but it is no longer necessarily an isomor- 
phism. The image of this inclusion is called the subspace of field-independent 
local horizontal forms. The quotient #£ ur (£) := tf" h ~f (£)/H^ p (M) is known 
as the space of conserved (local) p-currents. 

We also define the space of null local variational forms as the kernel 

of the projection fl*(F) —> f2*(£°°), which consists of all local variational forms 
that vanish on the PDE submanifold £°° . Therefore, we have the following 
exact sequence: 

&*<*(£) >n*'*(F) >Q*'*(£). (3) 

We denote the cohomology of the complex (f2*'°,dh) of null local horizontal 
forms by H* un (£). It should be noted that, provided the regularity conditions 
to be specified below are satisfied, null local forms are precisely those that are 
exact in terms of the Koszul-Tate part of the BV-BRST complex [35 ■ 

Given a vector bundle H — > M, we can analogously define the outer ends of 
the exact sequence 

t s (H)->r F (H)->r e (H), (4) 

where the space p£ (H) of null local sections of H consists of those that vanish 
on the PDE submanifold £°° and the space Ts(H) of on-shell local sections of 
H is the quotient. 

Now, let us consider the equation form e[4>] = as defining a local differential 
operator e: Tp{F) —¥ Tp(G) by the formula e[ip, <j)[ip]] = which corresponds 
to the smooth bundle map id x Pooe: J°°(F xf)-i J°°G. Consider the image 
Q°° = (id x p 00 e)(J°°(F x F)) C J°°(F x G). In this paper we are concerned 
with the case when Q°° is not necessarily all of J°°(F x G). In particular, we 
are interested in the cases when Q°° is itself a PDE in the sense given above. 
That is, there exists an involutive equation form z° : J°°(F x G) — > Z 1 such 
that Q°° is the zero set of irp x p x z° and a submanifold of J°°(F x G). The 
naturally associated local differential operator z° : Tp(G) — > Tp{Z 1 ) is variously 
known as a (stage-0) Noether operator, compatibility operator or redundancy 
operator, while we call Z 1 — > M the (stage- 1) Noether bundle. When there 
are no topological obstructions, the same construction can be iterated. Let 
Z° = G, (Z ) 00 = g°° and define (Z*) 00 = (tt f x p 00 z^ 1 )(J 00 (F x Z 1 - 1 )) with 
involutive equation form z % : J°°(F x Z l ) —> Z l+1 , with z 1 and Z l respectively 
the stage-i Noether operator and stage-i Noether bundle (which also naturally 
define local differential operators) . Provided there is no topological obstruction 
and the iteration terminate^] at i = r if (Z r )°° = J 00 Z r , the PDE system in 

1 Actually, at each stage, this procedure need only work on an open neighborhood of (Z r )°° , 
though then Z r+1 may need to be chosen as a non-linear smooth bundle. For simplicity, we 
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question is said to be stage-r irreducible. When r > 0, the PDE system is called 
a gauge theory, while when r = 1 it is said to be an irreducible gauge theory. The 
end point of this construction is a formally exact complex of local differential 
operators 

T F (F) r F (Z°) — Tf^Z 1 ) -^U • • ■T F (Z r ) ► 0, (5) 

where formal exactness |13[ 114) means that the following is an exact sequence 
of smooth bundles over M: 

Ps ° M2 '') ..J°°(fxZ r ) >Mx{*}. (6) 

To even have a hope of proving some stability results for the space of so- 
lutions of a PDE system, the PDE system itself must satisfy some regularity 
properties. Moreover, in the next section, we state a characterization of the 
characteristic cohomology groups in terms of the Noether complex, which only 
holds when appropriate regularity properties are satisfied. These properties 
are discussed in detail in [24j Sections 5.1, 6.4.2-3] and are summarized in the 
following subsections. 



2.2.1 Local regularity 

We say that a PDE system E°° c J°°F is locally regular if it is (i) stage- 
r irreducible and (ii) the corresponding differential operators e : J k : F G, 
z % : J ki (F x Z % ) — > Z l+X , where e and z % are respectively of orders k and can 
be chosen to be smooth bundle maps with Jacobians of constant rank and such 
that the Noether operators z l [V>,C] are linear in their second arguments, which 
we denote as z l [ip; £]. 

This requirement is essentially equivalent to that of [531 Section 5.1]. In 
particular, it allows us to conclude that any smooth bundle map / : J°°F — > F' 
that vanishes on the PDE manifold £ 00 (i.e., e[ip] — implies f[if>] = 0) must 
factor through the prolongued equation form: 

/: J°°F J-—> J°°(F x G) — ^— >■ F', (7) 

that is, f[ip] — f[ip;e[ip]}, where f[tp]^] is linear in the last argument. Further 
more, we can also conclude that any smooth bundle map g: J°°(F x G) — > G' 
that vanishes on the stage-0 Noether manifold Q°° = (Z ) 00 (i.e., g[ip, e[ip]] = 
for arbitrary ip) must factor through the prolongued stage-0 Noether operator: 

g: J°°(FxG)^^ x Z 1 ) — 9 -^G', (8) 

that is, g[ip, = g[tp; z°[ip; £]], where <7[^;C] is linear in the last argument. 
Similar remarks can be made about the higher stage Noether operators. In 
total, they are equivalent to the acyclicity of the Koszul-Tate complex in positive 
anti-field degree, which is used extensively in the BV-BRST analysis of gauge 
theories [24] , 

ignore these technicalities and work only with vector bundles over M. 
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2.2.2 Local linearizability 



In the question of linearization stability, we are interested in the relationship 
between the space of linearized solutions about a background exact solution 
ip: M F and a neighborhood of <p in the space of exact solutions. As such, 
the linearization of the PDE system in question must be well defined. 

Consider an arbitrary smooth 1-parameter family of sections fa : M — > F, 
with fa — (p + t-i/j + ^t, where ty t = 0(t 2 ). For convenience, we denote ipt = tijj + 
^t- We call a stage-r irreducible, locally regular PDE system locally linearizable 
about ip if (i) we can write, for fa as above and any 6 

e[fa] =e[y} + e v [i h ]- U[^t], (9) 

z i [fa;C} = z i X]-vU^;C% (10) 

where the local differential operators e v and z 1 [z^lC] = z % [(f;C]) are linear, 
while f<p[ij)t\ = 0(t 2 ) and yip[ipt]C] = 0(t), and (ii) the linear PDE system 
defined by e v [i/j] — is stage-r irreducible and locally regular, with its Noether 
complex given by 

r F (F) r F (z°) Tpiz 1 ) ■ ■ ■ v p(z r ) ► o. (11) 

We denote the PDE submanifolds defined by the linearized equation form e v [ip] = 
by 5™ C J°°F and £ v = ir k ^£™ C J k *F, where fc v is the order of e v . In- 
cidentally, linearization has made the local differential operators e v and 
field-independent. We can then naturally define the field-independent linearized 
Noether complex 

r(F) T{Z°) Y(Z X ) ■ ■ ■ T{Z r ) > 0, (12) 

which is a subcomplex of the field-dependent one above. 

The requirement of local linearizability also imposes a restriction on the 
allowed background sections ip: M — > F. We call such an admissible section ip 
a linearizable background. Denote by ru n (i ;l ) C T(F) the subset of linearizable 
background sections. Also, denote by S(£) C T(F) the subset of solutions, ip £ 
S(£ ) if e[<p] = 0. Finally, denote by Sy m (£) C S(£) the subset of solutions that 
consists of all linearizable backgrounds. All of these spaces can be topologized 
as subsets of T(F), with its natural smooth compact open (or Whitney) Frechet 
topology [25, 26 . We say that the PDE system in question is locally linearizable 
if the <Sii n (£) is an open subset of S(£). 

2.2.3 Local normality 

Finally, there is almost no hope of identifying linearized solutions with the full 
set of exact solutions if the non-linear system is of a higher order than the linear 
one. 

Let k, ki, k v and kf denote the respective orders of the differential operators 
e, z l , e v and z* . We say that a stage-r irreducible, locally regular, locally 
linearizable PDE system is locally normal if (i) the orders k — fc, = kf 
agree for all linearizable backgrounds ip £ S v (£) and (ii) the ranks of the linear 
bundle maps e v : J k F — > G and z % : J ki (F x Z % ) — > Z l+1 are respectively the 
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same as those of the Jacobians of the smooth bundle maps e: J k F — > G and 
z l : J k *{F x Z l ) -> Z i+1 . 

3 Deformations, cosymmetries, sources 
3.1 Consistent deformations 

Consider the linearized Noether complex [TT] Clearly, by linearity, the subspaces 
of null local sections, t £ (F) C T F (F) and t £ip (Z i ) C T F (Z i ) are preserved by 
the action of e v and zL . Therefore, we can define the null and on-shell Noether 
complexes as the top and bottom rows of the following commuting bicomplex: 

r £v (F) r Ev (z°) t £ , (zi) ■ • • r F (zn ► 

r F (F) — r F (z°) — ^ r F (z!) ■ ■ ■ T F {z r ) ► o, ( 13 ) 

(F) v £v (z°) r £v (z 1 ) ■ • • r F (zn > o 

where the columns are exact and the middle row is formally exact. We call 
the cohomologies of the rows of the above bicomplex, respectively, the null, 
(off-shell) and on-shell (stage-i) consistent deformations, denoted respectively 
by H\ et {e v ), H l Ae{ (F,e v ) and H l dci (e v ). Note that, since the Noether complex 
depends explicitly on the equation form e v , rather than just on the PDE sub- 
manifold £ v , we indicate the same dependence in the notation for consistent 
deformations. 

From the definition, an on-shell consistent deformation class [/] 6 if° ef (e v ) 
is represented by a local section / £ T F (Z ) = T(G). The local section / is said 
to be trivial if [/] = [0], which by local regularity means that it must be of the 
form 

m=e v [gm+h{i>;e v [i>}i (14) 

for some local section g £ T F (F) and some local differential operator h : T F (G) — > 
T F (G), with h[ip;Q linear in its second argument. 

Recall the linearization formulas |9] and [TU1 Let us define the leading order 
m > 1 and leading term of f v as 

UM = t m f^bp] + o(t m+1 ). (is) 

Similarly, let rtii > 1 and y$> m% ^ the leading order and leading term of . We 
also use the same notation for higher order terms in / and y i . Clearly, /^> m ^[?/>] 
and Vp™ 1 ^ [ip; C] are homogeneous in ip, of order m and mi respectively. Consider 
the following expansion of the exact Noether identity z [<fi; e [</>]] = 0, which holds 
for arbitrary sections <fi = f + tip 6 T(F): 

=>■ tz%[eM\ = t m z° v [f! p m) m+ty v [t^,eM}+O(t m+1 ). (16) 
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Completing the Taylor expansion and comparing the coefficients of powers of t, 
we find the condition z° [e^,^]] = at order 0{t), which was already a require- 
ment of local linearizability, and at order 0(t m ) that 

In other words, we find that represents an equivalence class [/p ] £ 

H^ ei (e v ) of on-shell consistent deformations. 

Of course, one could consider the cohomology spaces H^ e{ (M, e v ) of the field- 
independent linearized Noether complex f| 12[) . These field-independent consis- 
tent deformations are naturally included in those that are possibly field depen- 
dent, H l dc{ {M,e v ) C H^ ef (F 7 e v ). Of course, any non-trivial field-independent 
consistent deformation is still non-trivial on-shell, while any non-trivial null 
consistent deformation cannot be field-independent. Since in this paper we are 
chiefly concerned with non-linearities, all consistent deformations considered in 
the sequel will be field-dependent and of at least quadratic order. 

3.2 Cosymmetries and conserved currents 

For PDE systems that satisfy the regularity, linearizability and normality con- 
ditions given previously in Section 12.21 it is well known that there is a di- 
rect correspondence between the spaces HP nr (£ v ) of conserved local p-currents 
(which is the same as the space 

#cha~r P: \ £ ) I ' H dR P] '( M ) of on-shell closed local 
(n — p)-forms modulo field- independent ones) and the spaces of on-shell stage- 
(p — 1) cohomology of the complex formally adjoint to the linearized Noether 
complex 1111 We refer to the latter objects as stage-(p — 1) cosymmetries and 

denote their space by iJdosym(e v )) they are defined more precisely below. For 
simplicity, since we need only deal with the linearized Noether complex below, 
we restrict our discussion of this correspondence purely to the case of linear 
PDE systems. 

If h: J°°F — !> H is a linear differential operator, we define its formal adjoint 
h* as follows. For any vector bundle H — > M, denote by H* — > M its dual 
bundle and by H* = A n M ® H* the densitized dual bundle. There is a natural 
fiber- wise paring between sections r\: M — s> H and dual densities a* : M — > H*, 
n ■ a* : M -> M. The formal adjoint h* : J°°H* -> F* is defined as the unique 
linear differential operator that satisfies 

h[>p]-a* -ip-h*[a*]=dg h [iP,a*] (18) 

for arbitrary sections ip: M —> F and a* : M — > H* and some bilinear bid- 
ifferential operator G h [ip,a*}: J°°(F x H*) A rl_1 A/, where Q h is called a 
Green form associated to the adjoint pair h and h* . Note that we can consider 
g h as an F x iT-local horizontal form, Q h e Q, n - lfi (F x H*). Also, Q h is 
only defined up to the addition of an exact local horizontal form Oh[>P,(x*] ~ 
Qh\^,ot*\ + dfc%h/>, a*], where % is itself a bilinear bi-differential operator, so 
only the corresponding equivalence class [Qh] is well defined, though we may re- 
strict our attention only to representatives that are bilinear in their arguments. 
If h[tp; £] is a differential operator that is linear only in its second argument, we 
can still define its formal adjoint with respect to the second argument alone. It 
is again denoted by h*[ip; a*] and is also linear in its second argument, while an 
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associated Green form is denoted by Gh[i>>€,<x*] and is bilinear in its last two 
arguments. 

It is straight forward to verify that for a stagc-r irreducible, locally regu- 
lar, locally linearizable, locally normal PDE system, the formal adjoint of its 
linearized Noether complex is also formally exact. Also, since all the operators 
involved are linear, the subspaces of null local sections, Tg (F*) C T F (F*) and 
T £ip (Z r *) C T F (Z r *), are preserved. That is, we can construct the following 
commuting bicomplex, where the columns are exact, while the middle row is 
formally exact: 

^ (F* ) <-^- t £v (Z°* ) t £v (Z 1 *) ^— ■■•t £v {Z") < 

* 4- 0* ~ir 1* 

T F (F*) «-^— T F (Z *) Y F {Z 1 *) T F (Z r *) i 0. ( 19 ) 

F £v (F* ) +5- T £v (Z°* ) ^- T £v (Z 1 * ) <S- ■■■T £v {Z") i 

Finally, we dehne (for r > 0) null, (off-shell) and on- shell local (stage-r) cosym- 
metries to be elements of the kernel of the local differential operator z % * on, 
respectively, the top, middle and bottom rows of the above bi-complex; by 
convention we define z^ 1 = e v and Z^ 1 = F. That is, £ r * e T F (Z r *) repre- 
sents an on-shell local cosymmetry if z^~ r> *[ip\^' r *[^)\\ = whenever e v [xp] = 0. 
The cohomologies of these rows define equivalence classes of cosymmetries. A 
cosymmetry £ r * G T F (Z r *) is trivial if, respectively, 

(nuii) rm = z r v [t +u bp-,e v m, (20) 

(off-shell) f*M=*£[f +1 *M], (21) 
(on-shell) = z£[f + C*[4>; e^}}. (22) 

Two representatives are in the same equivalence class precisely when they differ 
by a trivial one. The spaces of equivalence classes of null, off-shell and on- 
shell stage-r local cosymmetries, respectively, by H^ osym (e v ) , H£ osym (F, e v ) and 
H^ osym (e v ). We call these equivalence classes rigid cosymmetries. The definition 
of a cosymmetry explicitly involves the equation form e v , rather than just the 
intrinsic PDE submanifold £ v , and the notation reflects that. 

Since all operators involved in the definition of the adjoint linearized Noether 
complex are field- independent, we can also define the field-independent adjoint 
linearized Noether complex 

r(F*) r(z *) r{z u ) r{z r *) < o, (23) 

whose cohomology spaces we denote by H£ osym (M, e v ). We have the natu- 
ral inclusion H^ osym (M,e v ) C H^ osyra (F, e v ) . Of course, any non-trivial field- 
independent local cosymmetry is still non-trivial on-shell, while any non-trivial 
null local cosymmetry cannot be field- independent. 

We conclude this section recalling the following known bijection between 
equivalence classes of higher stage cosymmetries and higher conserved currents. 
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We state the result only for the linearized equation, though a similar result 
works for non-linear equations as well. 

Proposition 1 (Generalized Noether's first theorem). For the linearized PDE 
system £ ip , the space of classes of conserved local currents is isomorphic to the 
space of classes of null local forms, which in turn is isomorphic to the space of 
classes of rigid on-shell local cosymmetries. That is, for < p < n, 

^cur {£<fl) = ^mlU^v) = ^cosym( e v)' (^4) 

Though the proof of Proposition [1] can be found in several places in the 
literature [Ml HUH- 

The name cosymmetry is meant to be evocative. For a variational PDE 
system (one defined by the Euler-Lagrange equations of a local Lagrangian) , the 
linearized equation form e v : Tp(F) Tp(F*) is formally self-adjoint, e* = e v . 
In that case, the adjoint linearized Nocther complex is identical to the formally 
exact complex of generators of (higher stage) gauge transformations (generalized 
Noether's second theorem). It is well known that the cohomology of higher stage 
gauge generators is identified with higher stage rigid symmetries, those of stage- 
being the ordinary symmetries. Stage-1 rigid symmetries have also been called 
reducibility parameters |12j . For non- variational system, there need not be any 
correspondence between gauge generators and Noether operators or between 
symmetries and cosymmetries, hence the need for the distinct terminology. 

In the case of stage-0 irreducible relativistic field theories and the Einstein 
or Yang-Mills gauge theories, rigid symmetries are known to correspond to so- 
called Killing vectors (or generalizations thereof 0[7]). By analogy, we can call 
the condition 

^(r-i)* K >* [v , ]] = and rm^z r v *li ir+1) *W (on-shell), (25) 
which defines a non-trivial cosymmetry, the (generalized) co-Killing condition. 

3.3 Null sources 

A null local p-source p is a local horizontal (n — p)-form, p G fl n ~ p '°(F), such 
that (i) it is horizontally exact, p = dhj for some local horizontal (n— p — l)-form 
j G £l n ~ p ~ 1 '°(F), and (ii) it vanishes on linearized solutions, that is, it pulls back 
to in n n -P>°(£ v ) along the inclusion £™ C J°°F, or simply p e {l n ^°(£ v ). 
A null local source p is said to be trivial if p — d^j where j itself vanishes on 
solutions. Two null local sources are said to be equivalent if they differ by a 
trivial one. We denote the space of null local p-source classes by HP rc (£ v ). 
The term source is meant to be evocative. Consider an equation of the form 

dj = p; (26) 

if the left-hand side can be considered as the divergence of a current, then the 
right-hand side should be considered the source (or source density) for that 
current, whence the exactness requirement. 

Null sources are clearly related to null local forms, as well as to characteristic 
cohomology. 
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Lemma 2. We have H™. c (£ v ) = H\ 



null 



(£ v ) = and for < p < n 



H^(£ v ) - H^(£ v ) S fl^- 1 ^)/^- 1 "^ W = ^C^)- (27) 



Proof. If follows directly from the above definition that for each null local p- 
source p, we can find a local horizontal (n — p — l)-form j such that d^j = p. 
Clearly j cannot be field-independent, unless it vanishes, and duj vanishes on- 
shell. But this means precisely that j is a conserved (p+ l)-current. Conversely, 
for each conserved (jp+ l)-current j, we can define p = d^j and easily check that 
p is a null p-source. Moreover, this correspondence respects equivalence classes. 
In other words, the classes of null p-sources are in bijection with the equivalence 
classes of conserved (p + l)-currents. The isomorphisms with H^fi^ip) and 

H"r*~ (£p)/-^dR P- C^O f°ll° w directly from definitions. 

The exactness of the variational bicomplex shows that any local horizontal 0- 
form is horizontally closed only if it is field-independent and constant. Therefore, 
any horizontally closed, null horizontal 0-form must be trivial. On the other 
hand, any null n-source must also be trivial, since the only exact local horizontal 
form in this degree is zero. In other words, H™ c (£ v ) — H®^^) = 0. □ 

Here is another convenient way to represent null sources. A local horizon- 
tal (n — p)-form p naturally defines a local differential operator p: T(M) — > 
Q, n ~ p (M). If p is a null local source, local regularity of the linearized PDE 
system implies that it can be written as p — p(e v ), where on the right-hand 
side we denote the local section e v £ Tp(G) is acted on by the linear local dif- 
ferential operator p: Tp(G) — > f2 n_p '°(F). We write this as p[ijj] = p[t/j',e v [ip]]. 
The local differential operator p[ijj; £] is not unique, since any other one of the 
form p[ip; -\-o~\tp; z°[Q] +r[ip, e^[*0], £], where a is linear in its second argument 
and r is bilinear and anti-symmetric in its last two arguments, would represent 
the same null p-source p[ip]. Local regularity implies that these possibilities 
exhaust the available ambiguity. We call a linear local differential operator 
p: T F (G) -> n*'°(F) trivial if it is of the form 



We have just observed that the space of equivalence classes of null p-sources 
Hi TC (£tp) is isomorphic to the space of equivalence classes [p] of linear local 
differential operators p: Tp(G) —> il n ~ p,0 (F) as above, modulo trivial ones. 

4 Linearization instability 

Consider a PDE system that is stage-r irreducible, locally regular, locally lin- 
earizable and locally normal, with the defining equation form and Noether com- 
plex Similarly, its linearized equation form and Noether complex about a 
linearizable background solution ip 6 iSi; n (£) are given by 1111 In this section, 
we shall refer to the space Su n (£) of linearizable background solutions simply as 
the (exact) solution space. We shall refer to the space S{£ v ) of solutions to the 
linearized PDE system at (p as the linearized solution space at (p. Both the exact 
and linearized solution spaces can be endowed with a topology as subspaces of 
the total space of sections T(F). The choice of topology on T(F) should be 
adapted to the problem at hand. The work on linearization stability that fol- 
lowed in the footsteps of [TB] identifies the space of solutions with the space of 



p[t, C] = C] + z°JQ] + t[V>; e v ty], Q. 



(28) 
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Cauchy data satisfying initial value constraints, which is endowed with a norm 
topology. Thus, the accompanying functional analytical steps that require com- 
pleteness introduce some weak solutions to the elliptic initial value constraints 
along with smooth classical ones. On the other hand, restricting our attention 
to a single Cauchy surface, rather than an open spacetime domain, does not take 
into account the possibility of singularity formation withing the domain, which 
excludes some initial data from the solution space. The seminal result on the 
global non- linear stability of Minkowski space [27J can be seen through the prism 
of linearization stability analysis. In that work, another normed topology was 
used on T(F). It would be interesting to also examine the same questions using 
a more natural Frechet (Whitney) topology [25, 26 in which T(F) is complete, 
without the need to introduce weak solutions, beyond the minimal attention 
this question has already received in the existing literature. 

We say that the background solution ip G Sy m (£) is linearization stable if 
(i) there is a neighborhood U C Sy sn {£) of ip that is homeomorphic to the 
linearized solution space SiS^) and (ii) each 1-parameter family of linearized 
solutions of the form tip, ip G S{£ v ) is mapped by this homeomorphism to a 
smooth 1-parameter family of exact solutions, tip i— ► <p t , where 4>t G <->lin(£) and 
<pt =tp + tip + 0{t 2 ). 

Obviously proving linearization stability is a difficult problem that must 
involve a significant amount of functional analysis and at present it has only 
been treated in detail for a few selected equations. On the other hand, it is much 
easier to show that linearization stability is obstructed. We devote the following 
subsections to explicitly exhibiting global geometric conditions on the manifold 
M and a background solution ip on it, which we call co- Killing conditions, that 
imply the existence of obstructions, also called linearization instabilities. 

4.1 Linearization obstructions 

Linearization stability is obstructed at a linearized solution ip G S(£ v ) if it 
cannot be extended to a smooth 1-parameter family of exact solutions of the 
form (p t = ip+tip + 0(t 2 ); if such a 1-parameter family does exist, then ip is called 
extendable. A stability obstruction is a function on the linearized solution space, 
Q: S(£ v ) — > V, where V is some vector space, such that Q(ip) = whenever 
ip is extendable. A linearization obstruction is said to be trivial if Q(ip) — 
for every linearized solution. The obstruction is said to be of order m if it is 
homogeneous of the same order, Q(tijj) = t m Q(ip). 

4.2 Deformation currents from null sources 

We are now ready to prove a theorem that relates null sources, consistent defor- 
mations and conserved currents. An easy consequence of it will be the existence 
of linearization obstructions, to be discussed in the next section. 

Theorem 3 (Deformation currents). Provided a linear PDE system £ v with 
equation form e v [ip] = is locally regular, there exists a natural bilinear mapping 
pairing a null local p-source class with a stage-0 consistent deformation class 
giving a conserved local p-current. 

j : Hl c {£ v ) x H° dcf (e v ) -> ^ur(^), (W. if}) [j P j]- (29) 
We call j p j the deformation current associated to p and f . 
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Proof. As discussed in the last two sections, provided the regularity condition 
is satisfied, each null local p-source class [p] can be represented by a linear local 
differential operator p: Tp(G) — > £l n ~ p '°(F), while a stage-0 on-shell consistent 
deformation class [/] can be represented by a local section / € Tp(G). We define 
the local horizontal (n — p)-form j p j £ il n ~ p '°(F) by the formula j p j = p(f); 
in other words, 

3 P jW\=p[tJW\l (30) 
for any section ip : M — s> F. It now remains to check that j p j is in fact conserved 
and that the map j is well defined on equivalence classes. 

We have already established that the horizontal differential of a null p-source 
representative p has the form 

dp[^C]=^;4[C]]+T[^e^[^] J C]- (31) 

The second term already vanishes on-shell because it is linear in the e v [ip] argu- 
ment. If we set £ = f[ip], the first term on the right-hand also vanishes on-shell, 
since a is linear in its second argument and ^[/[V 1 ]] vanishes on-shell by the 
defining property of a stage-0 on-shell consistent deformation. Therefore, j p j 
is in fact a conserved local p-current. 

If p is a trivial null local p-source, then it must take the form given in 
Equation ([251) . provided local regularity holds. We check that each of the three 
possible terms is a trivial conserved p-current after setting £ = f[ip]. The 
term dj[ip;(] is trivial because it is exact. The term a[ip; z® [f[ip]]] is trivial 
because z° [f[ip]] vanishes on-shell by the definition of a consistent deformation. 
Finally, the term r[ip] e<p[tp], f[ip]] vanishes on-shell because it is linear in the 
e,p[ijj] argument. 

If / is a trivial on-shell consistent deformation, then it must take the form 
given in Equation ([Til) , provided local regularity holds. We check that each 
possible term substituted for £ in p[ip; £] gives a trivial conserved local p-current. 
The second term gives p[ip; h[if>; e v [ip]]], which is linear in e^ip] and hence trivial. 
The remaining term p[ip; e v [g[ip]]\ is trivial for a slightly non-obvious reason: a 
deformation of the form e v [g [?/>]] necessarily comes from a local field redefinition 
of the form ip — > ip + tg[ip] + 0{t 2 ). Namely, recall that we can always write 
p[ip\ e^f?//]] = dfc["0] for some conserved local (p + l)-current k and consider the 
identity 

dfc[?/> + tg[ip]} = /# + tg[i/)];e v [il> + tg[^}}} (32) 
= tpty + t£#];e ¥ [£#]]] + p[if> + tg[il>];e v [il>]] (33) 
= e v [ 9 m + P W i^ e v W\]) + 0(t 2 ). (34) 

When expanded in powers of t, all coefficients on the left-hand side are exact. 
This shows that, up to the addition of the trivial conserved local p-current 
p^ [ip; e ¥ ,[V']], the term p[ip; e ¥ ,[g | [V']]] is exact and hence trivial. This concludes 
the proof. □ 

Remark 1. It is interesting to note that the deformation current mapping 

([p],[f])^\j P j], (35) 

defined in Theorem[31 identifies a selection criterion on consistent deformations. 
We could say that a consistent deformation / is conservative or p- conservative 
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when the deformation current j p j is trivial (on-shell exact). The linearization 
obstructions generated by deformation currents, as to be discussed in the next 
section, are necessarily trivial. 

4.3 Obstructions from deformation currents 

Consider the full non- linear PDE system £ defined by the equation form e[<fi] = 0, 
which we take to be locally regular, linearizable and normal. Also, consider a 
linearizable background solution cp 6 5ii n (£). Recall that, a smooth 1-parameter 
family of solutions, e[<fit] = 0, of the form cf> t — ip + ipt satisfies Equation ([5]), 
that is, 

e v [i>t] = UbPt]- (36) 

If ipt — tip + 0(t 2 ), then ip is a linearized solution, e v [ip] = 0, that can be 
extended to a smooth 1-parameter family of exact solutions, for instance (pf 
Recall also that the non-linearity f v defines a representative of a consistent 
deformation class [/<»] e H^(e v ) given by the Taylor expansion 

UM=t m f {m) m+0(t m+1 ) (37) 

for some m > 1, where f( m ^[ip] is homogeneous of degree to in -0. The main 
result of this paper, to be proven below, is that each deformation current defined 
by Z^'™- 1 canonically defines a linearization obstruction. 

As an intermediate step before formulating the main result, we need to 
introduce a grading by degree of polynomial dependence on field bundle sections 
on the spaces of local sections and local differential operators. Let ^Tp(H) C 
Tp(H) denote the subspace of local sections of a vector bundle H — > M of 
homogeneous polynomial degree I, that is, such that h G ^Tp(H) only if h[tp] 
depends polynomially on i\> and its derivatives, and is homogeneous in ip of 
degree I. Note that the horizontal differential dh, as well as all the differential 
operators in the linearized Noether complex 111! the adjoint linearized Nocthcr 
complex [T51 and the complex Tp(H) — > Tp(H) — > T£(H) all preserve subspaces 
of homogeneous polynomial degrees. This means that any of the cohomology 
spaces that we have defined also have subspaces of homogeneous polynomial 
degrees; the corresponding subspaces are denoted by ^H^(-) C H*(—). The 
representatives of each of these subspaces can always be chosen of the same 
homogeneous polynomial degree. 

We are finally ready to formulate and prove the main theorem of this paper, 
which relates the de Rham cohomology H^ R (M) and the space H* osym (£ v ) of 
rigid on-shell cosymmetries with linearization obstructions at if, and thus with 
potential linearization instability. 

Theorem 4. If ip € S(£ v ) is a linearized solution that can be extended to a 
smooth 1-parameter family of exact solutions, then it must necessarily satisfy 
the conditions Q l p (ip) = 0, < p < n, < I, where 

Q l p ■ S(€ v ) -> ^Hl mym {e v y g) K R P (M) (38) 

are linearization obstructions of order I + m defined by the non-linearity f v , 
where * denotes the linear dual. 
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Proof. First, recall the isomorphisms (0 < p < n) 

^cosym( e v) = ^cur (^¥>) = ^src(^v) (39) 

and note that they preserve subspaces of homogeneous polynomial degree. That 
is, we can consider a null local p-source class [p] G ( l ' Hg IC (£ v ) , equally well, to be 
an element [p] G ^^cosym( e v)- Using the on-shell consistent deformation class 

Lfp"^] G ^dcf( e v) ano - the ma P j defined in Theorem [3J we obtain a conserved 
local p-current class 

km, i4 m) ]) = y pj ^} e (i+m) ^ ur (^) = a+m) ^(^)- (40) 

The representative j .<,„> G il™~ p '°(F) is a local horizontal (n — p)-form and 
hence naturally defines a differential operator j ,< m ) : — > Vl n ~ p (M). By 
construction, we know that j .(,„> hp] is of homogeneous polynomial degree l+m 
and that dj ,„„ [•</-] = if ■(/' is a linearized solution. Hence it represents the 
de Rham cohomology class in [j p y(™> ["0] ] G H^ P (M). Finally, we define the 
map Q l p by the formula 

[p]-Q l M) = [3 pJ ^W\l (4i) 

where the dot on the left-hand side stands for the natural pairing between 
( °ffc P osym(e v ) and its dual «ff? osym (e^)*. 

Since Q p (ip) is built out of a bilinear pairing between the differential operator 
p and ftp hp], the order Q l p (ip) must be exactly / + m. 

To conclude the proof, it is sufficient to show that [p] ■ Q l p {ip) = [0] for ip that 
is extendable to a smooth 1-parameter family of exact solutions cpt = <fl + ipt> 
with ip t = tip + 0(t 2 ). This family of exact solutions will satisfy the equation 
e <pbPt] = ftpH>t]- Recall also that we can find a conserved local (p + l)-current 
k G O n -P-h°(F) such that p[^,e v [^]] = dfc[*] for arbitrary * G T(F) and, 
in particular, this formula is still valid if we set ^ = ip t . We then have the 
following equalities in terms of de Rham cohomology classes: 

[p^fJAA]}} =t (l+ ^[p[^f^m}+0(t l + m+1 ) (42) 
= [phfc;e^]]] = [d*#t]] (43) 
= [0]. (44) 

Note that comparing the coefficient of t l+m in these equations gives precisely 
the desired equality [p] ■ Q l p {ip) = 0. □ 

Remark 2. The above theorem only shows how to canonically construct po- 
tential linearization obstructions. It does not mean that the obstruction Q l p is 
necessarily non-trivial. For instance, if one can show that the conserved current 
j ,(m) is trivial, [j f < m >] = [0] <G HP m (£ v ) (equivalcntly, that the consistent 

deformation /( m ) is p-conservative, as defined in the preceding section), then 
[p] ■ Q l p {4>) — [0] is trivial and does not pose any obstruction to extending lin- 
earized solutions to exact ones. Further work must be done in on a case by case 
basis to show that the obstruction Q l p is non-trivial. Alternatively, one can use 
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the triviality of [p] ■ Q l p (ip), checked in the way just given above, as a condi- 
tion to select non-linear consistent deformations like that will not create 
linearization obstructions (at least not of the kind constructed in Theorem [4]). 

Theorem 2] is a significant generalization of the results obtained in the lit- 
erature that followed up the initial work of |28l 116) . The known linearization 
obstructions have only been obtained in the case of either compact manifolds 
for equations in Riemannian geometry or in the case of manifolds with compact 
Cauchy surfaces for relativistic gauge theories. In light of our result, these sit- 
uations are immediately recognizable as obstructions coming respectively from 
H n (M) ^ and from H n ~ 1 (M) ^ 0, and also why only gauge theories are sus- 
ceptible in the latter case. Since the equations considered in the literature have 
been stage-0 or stage- 1 irreducible, our result also shows why similar obstruc- 
tions are absent for manifolds with H n ~ p (M) ^ andp > 1. To see obstructions 
due to these lower cohomology groups, one needs to consider non-linear equa- 
tions with higher stage reducibility and they are relatively infrequent. Such 
gauge theories do appear in some specialized physics literature, with non-linear 
p-form electromagnetism and some supergravities as examples. They have also 
recently attracted significant attention in the mathematics literature as higher 
gauge theories ISO] • 

Finally, the calculations involved in identifying these linearization obstruc- 
tions have been rather cumbersome. This is especially the case for relativistic 
gauge theories, where the necessary calculations lose much of their geometric 
character due to non-covariant restriction to some initial data surface. Our re- 
sult, on the other hand, was obtained completely geometrically and explains the 
covariance of the final result. Famously, for Einstein equations, linearization 
obstructions are only present if the background solution possesses non-trivial 
Killing vectors (or satisfies the Killing condition). In particular, in the existing 
work on relativistic gauge theories, significant effort was needed in each case 
to obtain the analog of the Killing condition. Our result, on the other hand, 
shows that this condition is precisely Equation [55] defining a non-trivial cosym- 
metry, which for obvious reasons we have also called the (generalized) co-Killing 
condition. 

The relation of our result to previous work is discussed in more detail next. 

4.4 Integrated charges, relation with previous work 

Note that the linearization obstruction [p] ■ Q p (ip) constructed in the preceding 
section is not strictly of the kind that has been constructed in the existing liter- 
ature on Einstein, Yang-Mills and related equations. While both formulations 
pass through conserved currents, we have given the obstructions as valued in the 
de Rham cohomology H^ R (M) classes of these currents, while the usual formu- 
lation gives them in terms of corresponding integrated charges. An integrated 
charge is obtained by integrating a conserved current over a closed compact 
submanifold. 

Our formulation can give rise to integrated charges as well. If E C M is a 
closed compact submanifold, then 




(45) 
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is the corresponding integrated charge. The charges obtained in this way are 
not all independent. Each charge depends only on the (singular) homology class 
[£] G i?*(M). Moreover, linear combinations of homology classes lead to linear 
combinations of charges. Therefore, to obtain a set of independent charges, 
we need to pick a basis [Si] for H*(M) and pair each basis element with [p] ■ 
Q p (ip). But this amounts to no more than composing Q l p with the isomorphism 
iJ* R (M) = 0,M bl defined by the basis [£»], where the b l = dim H l (M) are 
the Betti numbers of M . The fact that this map is an isomorphism is simply 
a restatement of Poincare duality |31) . So, an integrated charge over a non- 
trivial closed compact submanifold £ C M is simply a witness to the existence 
of a non-trivial cohomology class in H*(M). In other words, our formulation of 
linearization obstructions is equivalent to the usual one. 

4.5 Asymptotic boundary conditions 

On the other hand, our formulation is more convenient in the discussion of 
asymptotic boundary conditions on non-compact manifolds and corresponding 
non-standard de Rham cohomology. For instance, if M is non-compact, we may 
consider a suitable subspace of T(F) C T(F) of the space sections of the field 
bundle F — > M, along with a corresponding refinement fi fc (M) C f2 fe (M) of 
the de Rham complex, with d£l k (M) C fr +1 (M). These subspaces could be 
selected by imposing some asymptotic boundary conditions at the open ends 
of M. An extreme example would be require all sections and forms to have 
compact support. The cohomology H^ R (M) = 7J(fi*(M),d) could be different 
from the standard H^ R (M). Consider, with respect to the linearized PDE 
system, a conserved local current class [k] and the corresponding null local 
source class [p] = [d^k]. If the local horizontal forms p and k can be chosen such 
that the boundary conditions on ip G f (F) imply that p[ip], k[ip] e fl*(M). The 
construction of the linearization obstruction {p\-Q p {\p) still works, but the result 

is now valued in the non-standard de Rham cohomology H*j. R (M) rather than 
H%j,(M). As before, it is still a non-trivial problem to check that the resulting 
potential linearization instability Q l p is non-trivial. However, if the boundary 
conditions are chosen such that H*(M) = 0, then the linearization obstruction 
yielded by this construction is necessarily trivial. 

The above logic appears to be the reason behind the absence of linearization 
obstructions for the common choice of asymptotically flat boundary conditions 
for the Einstein equations [3j 0] . 

5 Examples 

The PDE systems studied in the physics literature are mostly variational (com- 
ing from classical Lagrangian field theories) . These are also the kind of systems 
analyzed in previous work on linearization instabilities. Our analysis is appli- 
cable to a more general class of PDE systems, including non-variational one. 
It would be nice to have explicit examples of theories from each of the corners 
missed by the existing literature. On the other hand, it is rather easy to pro- 
vide examples of non-variational PDE systems by taking a variational system, 
e[if>] — 0, and pre-composing with an arbitrary differential operator, e[g[r/]] = 0. 
The resulting system will generically no longer be variational, but will possess 
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essentially the same Noether complex as the original one. Thus, the essen- 
tially new examples that we give below are restricted to 1-dimensional systems 
(ODEs) and higher gauge theories. 

5.1 Stage-0 irreducible systems 

5.1.1 ODEs 

Consider an ordinary differential equation (ODE) that is defined on a circle, 
M = S 1 , and scalar valued, G = F = M xR: 

= iM (46) 

with <f) = dcp/dt and t a coordinate on M . This is essentially the 1-dimensional 
particle equation with cyclic time and force term /, which, say, is homogeneous 
of polynomial degree > 1. If we linearize about ip — 0, the linearized equation 
is just the free particle equation (f> — with cyclic time. Its only solutions are 
constants, S(£ v ) — {ip(t) \ ip(t) = tp } — K 1 - It is straightforward to check that 

= dt is a non-trivial cosymmctry with corresponding conserved 1-current 
k[ip] = ip (the momentum) and null source p[ip] = ipdt. 

If the non-linearity is f v [ip] — f[ip] — ip , the corresponding deformation 
0-current is 

j p ,^[ip}=ip 2 dt, (47) 
which is easily seen to be non-trivial. In fact, the integrated charge 

{[S 1 ], Q p m - / 3 P ,vW\ = ctP 2 o 0, (48) 

with c the circumference of S , for any linearized solution ip(t) — ipo 0. 
Clearly, the zero set Q p (tp) = consists of a single point, {ip{t) \ ip(t) — 0} = MP. 

On the other hand, if we consider £* = ipdt, j[tp] = ^tp (the energy) and 
p[tp] = tpip dt, the corresponding deformation current is trivial: 

Jp ,^]=^ 2 dt=idVA (49) 
In other words, the force term /[</>] = 4> 2 is energy-conservative. 

5.1.2 Semilinear elliptic PDEs 

A very similar situation occurs on higher dimensional compact manifolds, say 
M — S 2 , with Laplace-type semilinear equations, F = M x R, G = F* = A 2 M, 

*(A0 + l(l + 1)0) =f[ct>}, (50) 

where we interpret S 2 as the standard, round unit sphere, * is the Hodge star, 
A = *d*d + d*d* is the Laplacian on it and / > is an integer. The lin- 
earized equation about ip — is *(A^ + 1(1 + = and its solution space 
is well known to be of dimension 21 + 1 and spanned by spherical harmonics, 
S(£ v ) = M. 2l+1 . A non-trivial cosymmetry is given by a spherical harmonic 
£* = ^m- The corresponding conserved 1-current is k m [4>] = *(Yj n d(j) — (pdY^) 
and corresponding null 0-source is pm[4>] = *^m(^0 + + ^)4>)- 
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If we set ftp[ip] — f[ip] — tp , we get the non-trivial deformation current 

j Pm! ^ = *Y^ 2 , (51) 

with integrated charge 

l\S\Q m m= I *Y^ 2 . (52) 
Js 2 

In other words, in order to satisfy all linearization obstructions of the form 
Qmiip) = 0, — i < m < /, the decomposition of ip 2 into spherical harmonics 
must have vanishing coefficients. These conditions are only satisfied by the zero 
solution, these obstructions drop the dimension of the linearized solution space 
from 21 + 1 to 0, {tp | tp = 0} ^ M°. 

5.1.3 Other examples 

Other examples considered in the existing literature concern some problems from 
Riemannian geometry on compact manifolds, such as the equation describing 
metrics of constant scalar curvature (32]. The identified obstructions can also 
be computed using our main result. 

5.2 Stage- 1 irreducible systems 
5.2.1 Einstein equations 

Einstein's equations describe the dynamics of the gravitational field, a Lorentzian 
metric, on a manifold M either in vacuum or, with appropriate modification, 
in the presence of matter. A Lorentzian metric g 6 T(S 2 T* M) is a section of 
the bundle F = S 2 T*M — >■ M of rank-2 covariant symmetric tensors. Detailed 
formulas necessary to exhibit the structure of the equations, the non-linearity 
and the linearization obstructions are rather lengthy and, for our purposes, not 
particularly illuminating. The linearization instabilities of Einstein equations 
have been studied extensively, they in fact gave birth to this subject, and all 
the relevant details are summarized in the introductory sections of [lj. Below, 
we merely indicate how they fit into the framework of our main result. 

The linearized Nocther complex for Einstein equations, as well as its formal 
adjoint, about a background metric tp = g, are given by 

T F {S 2 T* M) y Tp(T*M) — ^— > 0, (53) 

T F (S 2 T*M) < Tp(T*M) <r^- 0. (54) 

Note that we have identified F* = F, and the same for the other tensor bundles, 
using g to raise and lower tensor indices as well as to construct a canonical 
volume density. Here, L is the differential operator of the linearized Einstein 
equations, also known as the Lichnerowicz operator. The operators B and K 
correspond, respectively, to the linearized Bianchi identity and Killing equation. 
In local coordinates, they are 

{B[h])i = V j hij, (55) 
(KMh = \(y i v J +V J v i ), (56) 
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where V is the Levi-Civita connection with respect to g, which is also used 
to raise and lower tensor indices. Solutions to the Killing equation are called 
Killing vectors; in our framework they are the rigid stage- 1 cosymmetries of the 
linearized Einstein equations. The conserved 2-currents corresponding to these 
Killing vectors are known as Abbott-Deser fluxes [33 and the corresponding 
null 1-sources do not have a name in the literature. On the other hand, the 
corresponding deformation 1-currents, constructed using the quadratic term in 
Einstein's equations, are known as the Taub conserved currents 34 . The con- 
nection between the presence of non-trivial Killing vectors and the vanishing 
of their Taub charges as a linearization obstruction was first noticed by Mon- 
crief [6]. 

5.2.2 Yang-Mills equations 

The notational background for this section is given in Appendix [X] In Yang- 
Mills theory [33] , the basic dynamical field is a semi-simple Lie algebra g- valued 
1-form a G SI 1 (M,g), so F — g <S> A 1 M. Where M is a manifold of dimM = n 
and endowed with a (pseudo-)Riemannian metric. It is a stage-1 irreducible, 
non-linear deformation of the linear Maxwell theory. Its Lagrangian density 
with the leading order non-linearity, as a local variational form, is 

C = -^{d h a A *d h a) + i(d h a A *[a A a}) + 0(a 4 ). (57) 

The linearized equations, about (p = a = 0, and the leading order non-linear 
consistent deformation are obtained by a vertical variation of the Lagrangian 
density: 

d v £ = — i(d v dha A *dua) + i(d v dha A *[a A a]) + (d^a A *[d v a A a]) (58) 
+ 0(a 4 ) 

= — -(d v Q A *5hdha) + ^(d v a A *<5h[a A a]} + (d v a A [a A *dh«]} (59) 
+ d h (-.-) + 0(a 4 ). 
We can read off the direct and adjoint linearized Noether complexes as 

n i,o (F) fl) e ^** hdh , &'°(F, a) + ^°(F, a) ► 0, (60) 

n 3 '°(F, g) ( C ^* Ahdh g) Q >°(F, g) < 0, (61) 

while the leading order consistent deformation is 

fW[a] = f[a] = *6-[a A a] + [a A *da]. (62) 

A stage-1 rigid cosymmetry is any g- valued 0-form e G £l 0,a (F, g) such that de = 
0. The corresponding conserved 2-current and null 1-source are, respectively, 
k = (e A *da) and p — (e A *8da) (up to signs). The resulting deformation 
current, which we can check is not off-shell closed, is 

j pJ [a] = (e A [a A *da]>, (63) 
dj p j[a] = -{e A [a A (d*da)]), (64) 
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where the neglected term is trivial because of the identity d*5 = 0. 

The above construction gives an obstruction Q e for each rigid cosymmetry 
e e n 0,0 (F, g) valued in (M). These are precisely the obstructions that 

were previously obtained by Moncrief [7], where he checked that they are non- 
trivial and also sufficient. Moncrief, like in other existing literature, implicitly 
used the existence of a compact Cauchy surface as a witness to the non-triviality 
o{H^(M). Obstructions appear also at non-vanishing backgrounds tp = A, 
where the rigid cosymmetries are geometrically identified with A-parallel g- 
valued scalars. If A is interpreted as a connection on a principal bundle, these 
parallel scalars are intimately connected with its the holonomy group. 

5.2.3 Chern-Simons 

The notational background for this section is given in Appendix [21 In Chern- 
Simons theory [36], where dimM = 3, the basic dynamical field is also a semi- 
simple Lie algebra g- valued 1-form a € ri 1 (M,g), so F = g(g>A 1 M. It is a stage- 
1 irreducible, non-linear theory with Lagrangian density, as a local variational 
form, 

£= i(aAd h a) - -(a A [a A a}). (65) 

The linearized equations, about <p = a = 0, and the non-linear consistent de- 
formation are obtained by a vertical variation of the Lagrangian density: 

d v £ = (d v a A d^a) - (d v a A [a A a]) + d h (- • •)• (66) 

We can read off the direct and adjoint linearized Noether complexes as 

fli.0(F jfl ) e ^ dh ) n 2 -°(F,fl)-^n 3 -°(F,g) >0, (67) 

r! 2 '°(F, ) ( &>dh — n 1 ' ^, g) n°>°(F, ) < o, (68) 

while the leading order consistent deformation is 

4 2) [a]= /[«] = [<* A a]. (69) 

Just as in the case of Yang-Mills theory, a rigid stage- 1 cosymmetry is any q- 
valued 0-form e € fl° : °(F, g) such that de = 0. The corresponding conserved 
2-current and null 1-source are, respectively, k = (s A a) and p = (e A da) . The 
resulting deformation current, which we can check is not off-shell closed, is 

j pJ [a] = (eA [a A a}), (70) 
6j Pjf [a] = -{eA[aA(da)]>. (71) 

The geometric interpretation of the rigid cosymmetries, and their connection 
with the holonomy group, remains the same as in Yang-Mills theory. 

The solutions of the Chern-Simons equations constitute flat connections on 
the corresponding principal bundle over M (in the above simplified setting, 
the principal bundle is trivial). The moduli space of flat connections, that is, 
the solution space we have been considering modulo the gauge transformations 
a — » a + duj, has been studied extensively. It is known that this moduli space 
is not a smooth manifold, because it may have quadratic algebraic singularities 
at connections with non-trivial holonomy groups |37j . That result is entirely 
consistent with the obstructions computed above. 
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5.3 Other examples 

Other irreducible gauge theories that have been considered in the existing liter- 
ature include Einstein- Yang-Mills [5] and classical N = 1 supergravity [5]. The 
obstructions and (co-)Killing conditions identified in these theories are repro- 
duced by our framework and have been checked to be non-trivial and sufficient. 
See [3 [5] for a historical discussion and detailed references. 

5.4 Stage > 1 irreducible systems 

5.4.1 Non-abelian Freedman-Townsend 2-form 

The notational background for this section is given in Appendix [21 General- 
izations of Maxwell equations to p- forms (abelian p- forms) provide a standard 
source of examples of higher stage irreducible linear PDE systems. Their con- 
sistent deformations (non-abelian p- forms) provide examples of corresponding 
higher stage irreducible non-linear PDE systems. 

On a manifold M of dim M — 4 that is endowed with a (pseudo-)Riemannian 
metric, the only consistent, non-linear deformation of abelian 2-form field the- 
ory is the Freedman-Townsend model |38j . The basic dynamical field is a semi- 
simple Lie algebra g-valued 2-form /3 e Q 2 (M, g), so F — g ® A 2 M. Its La- 
grangian density, with the leading order non-linearity, as a local variational 
form, is 

C\p\ = -i(d h /3 A *d h /3) - A [*d h /3 A *d h /3]) + 0(/3 4 ). (72) 
o 4 

The linearized equations, about tp = /3 = 0, and the non-linear consistent defor- 
mation are obtained by a vertical variation of the Lagrangian density: 

d v £[/3] = - J(d v d h /3 A *d h /3) - J(d v /3 A [*d h( S A *d h /3]} (73) 
+ i(/3A[*d h /3A*d v d h /3])+0(/3 4 ) 
= -|(d v /3 A *5 h d h f3) - i(d v /3 A [*d h /3 A *d h /3]) (74) 
+ i(d v /?A*<5 h [/3A*d h /?])+0(/3 4 ). 
We can read off the direct and adjoint linearized Noether complexes as 

n 2,0 = > tffi_±+ S p,0_* L +tffi (75) 

while the leading order consistent deformation is 

4 2) [/?] = /[/3] = [*d/3 A *d/3] - *5(2[*d/3 A /?]). (77) 

A rigid stage-2 cosymmetry is any g- valued 0-form e G Sl°'°(F, g) such that de = 
0. The corresponding conserved 3-current and null 2-source are, respectively 
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k = (e A *d(3) (a 1-form) and p = (e A *Sd(3) (2-form) (up to signs). The 
resulting deformation current, which we can check is not off-shell closed, is 

j P ,M = (eh[*dpA*dp]), (78) 
dj P ,M = <£A[d*d/?A*d/3]>, (79) 

where the neglected term is trivial because of the identity d*5 = 0. 

Therefore, for each such rigid stage-1 cosymmetry e, we have found a second 
order linearization obstruction for the background solution ip = /3 = of the 
Freedman-Townsend model on any (pseudo-)Riemannian 4-manifold with non- 
vanishing iJj R (M): 

[e]-Q p :S(£ v )^Hl K (M). (80) 

An example of a Lorentzian, globally hyperbolic spacetime is the Schwarzschild 
black hole spacetime, whose topology is M 2 x S 2 . This observation appears to 
be new. We have not explicitly checked that these obstructions are non-trivial, 
but that is likely to be the case, by analogy with the Yang-Mills and Chern- 
Simons examples. It would be interesting to obtain an interpretation for these 
rigid cosymmetries in terms of the higher bundle interpretation of higher gauge 
theories |2"9"] . 

5.4.2 Other examples 

Higher gauge theories have attracted a lot of attention recently in some mathe- 
matical literature (for instance, sec [29, 30, 39 ). Many of these theories involve 
Lie algebra g- valued p-forms as dynamical fields and exhibit non-linearities sim- 
ilar to those of Yang-Mills, Chern-Simons and Freedman-Townsend theories. 
In addition, N > 1 classical supergravity theories also involve dynamical p- 
forms, non-linearly coupled to other fields .40!. It appears that the question of 
linearization stability or instability of particular background solutions in these 
theories has yet to attract any attention. 

6 Discussion 

We have considered the question of linearization instability for a large class of 
non-linear PDE systems, which includes relativistic Lagrangian field theories, 
with both irreducible and reducible gauge theories among them. This class 
consists of all (not necessarily Lagrangian) PDE systems, for which a Noether 
complex can be defined. 

The question of linearization stability has two aspects. One is to identify ob- 
structions, which prevent an arbitrary linearized solution from being extended 
to a family of exact ones. The other, once sufficiently many obstructions have 
been identified, is to show that no further obstructions exist. We have concen- 
trated only on the first aspect. Moreover, we have restricted our attention only 
those obstructions valued in the domain manifold's de Rham cohomology, with 
the cohomology representatives being constructed locally out of the linearized 
solutions. Though, within this class of obstructions, we likely give an exhaustive 
classification. (It is a matter of connecting the result of Theorem [3] with the 
consistent deformations coming from the non-linearity at leading, as was done 
in Theorem SI or higher orders). Interestingly enough, it captures essentially all 
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linearization obstructions that have been identified in the work following up the 
original articles Brill & Deser [35] and Fischer & Marsden [IB]. In many cases 
these known obstructions have also been shown to be sufficient [TJ [5] , at least 
when concentrating on initial data and ignoring blow-up singularities. 

Our main result (Theorem [4j provides a streamlined, unified way of iden- 
tifying such obstructions. Also, both the results and intermediate calculations 
are everywhere geometric and covariant with respect to the PDE domain. This 
is an improvement over previous calculations, which treated individually each 
PDE system of interest and often required some non-geometric intermediate 
steps. 

The conditions identifying linearization unstable backgrounds in known ex- 
amples have all been, in a sense, generalizations of the Killing condition (exis- 
tence of non-trivial Killing vectors) in general relativity. Our main result puts 
this observation into a more general context. We term the most general form of 
this condition the co-Killing condition (with the prefix co- implying that it is 
defined by the Noether complex, rather than the complex of gauge generators, 
which are identical in Lagrangian theories). This condition corresponds to the 
presence of non-trivial rigid (higher stage) cosymmetries (Section 13.21) . In La- 
grangian theories (by a generalization of Noether's second theorem), these are 
dual to rigid (higher stage) symmetries, of which Killing vectors are in fact an 
example. 

Our result also clarifies the role of the non-trivial topology of the PDE do- 
main. The known linearization obstruction have so far appeared mostly on 
PDE domains (say of dimension n) that are compact (H^ R ^ 0) or with a 
compact initial data surface (-ff^R ^ 0)- No particularly clear statement has 
been made about non-compact domains. Now, from Theorem |4j we know that 
any non-trivial de Rham cohomology class in -ff^R can generate a linearization 
obstruction. However, the lower degree cohomology classes come into play only 
in the presence of rigid higher stage cosymmetries. In particular, since reducible 
gauge theories have not been analyzed in the literature on linearization instabil- 
ities (Einstein, Yang-Mills and related theories are irreducible gauge theories), 
this explains why only the -H^r ancl ^dR cohomology classes have played a 
role. Furthermore, our result shows that linearization obstructions can appear 
even in the absence of non-trivial topology, as long as the de Rham complex 
develops non-trivial cohomology when restricted by some asymptotic boundary 
conditions on open manifolds (Section I4.5[) . 

Further, the intermediate result of Theorem [3] may be of interest in its own 
right. It shows that consistent deformations of linear PDE systems relate on- 
shell conserved currents of different degrees. Namely, a consistent deformation 
and a conserved (p + l)-current bilinear ly define a p-current, which we have 
called a deformation current. We have named conservative those consistent 
deformations that have trivial deformation currents. It may be interesting to 
analyze the subset of conservative deformations and its relation to the general 
problem of adding non-linear interactions to linear field theories |41) . In par- 
ticular, it is interesting to investigate whether the rank (in either argument) of 
the deformation current mapping is expressible is some known invariant of the 
linear PDE system. 

Finally, in the problem of deformation quantization of classical theories, 
the differential geometry of the phase space plays a crucial role. When the 
phase space is a smooth manifold, the Fedesov construction essentially solves 
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the problem [32]. In classical field theories, the phase space can be identified with 
the space of solutions. Non-smooth, singular points (precisely those solutions 
that are linearization unstable) then pose an obstacle to quantization. The 
situation with deformation quantization in the presence of such singularities is 
much less clear, though some work has been done in that direction [43l l44l 145] , 
We hope that the results of this paper could serve as tools in this program. 
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A Dynamical Forms 

Below, we collect useful formulas for explicit calculation with theories where the 
dynamical fields are differential fc-forms on the manifold M of dimM = n. We 
take the field bundle to be F = A k M, often the equation bundle and the Noether 
bundles are bundles of forms over M as well, say G = K l M and Z % = A ki M. The 
space of F-local section can then be identified with the space of local horizontal 
forms of the same degree, T F (F) = O fc <°(F), T F (G) ^ ft l '°(F) and I>(Z 4 ) = 
n k ^°(F). Dual densities can be identified with forms of complementary degree, 
F* = A"~ fc , where the natural fiber- wise pairing is realized as 

A-^ = AA/i, (81) 

with A G r(F) and fx G T(F*). The de Rham differential d extended to act on 
local sections is none other than the horizontal differential dh from the varia- 
tional bicomplex. Using the algebra of differential forms, it is easy to find its 
formal adjoint: 

dA A /i — A A (— )' A ' +1 d/i = d(A A fi), (82) 
which means that the formal adjoint of d h : fl lfl (F) -> Q, l+1 >°(F) is dJJ = 

(-) /+1 d h : n n - l - lfi (F) -> n n - k . 

When the manifold M is endowed with a (pseudo-)Riemannian metric g, 
we can define the corresponding Hodge dual operator *: il l '°(F) — > fl n ~ l '°(F). 
Recall the following identities [46] : 

Ai A *A 2 = A 2 A *Ai, (83) 

**A = £| A |77, With 6 k = (_)fc(n-fc) + (n- S )/2 ) (g4) 

where s is the signature of g. If n = 4 and the metric has Lorentzian signature 
( —+++), s = 2, then = (— ) k+1 . With Hodge duality available, it is also 
possible to identify the dual densities of /-forms with /-forms themselves with 
the Hodge fiber-wise pairing 

A-^ = AA*^. (85) 

The formal adjoint of dh with respect to the above paring is the horizontal 
dc Rham codifferential denoted by <5h (and by 5 when not acting on local hori- 
zontal forms): 

dhA A*/i-/iA *<5h^i = dh(^i A w). (86) 
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It satisfies the identity 5hA = (— )' A 'e| A |dh, so that 6% = and dh*#h = <5h*dh = 0. 

For calculations that include vertical forms, it is convenient to extend the 
Hodge dual to all local variational forms so that 

*(u A A) = v A *A, (87) 

where v is any purely vertical form and A is any purely horizontal one. This way, 
the Hodge dual commutes with the vertical differential, d v * = *d v . Moreover, 
the Hodge dual paring satisfies the identity 

(wAA)A *(vr A /i) = (-) |A|M w A tt A (A A (88) 
_ (-JCW+IwDM^aw A(/xA*A) (89) 
_ (^lAIM+MM+MM^ A M ) A A A ) ; (go) 

whenever w and tt are purely vertical, while A and fi are purely horizontal. 

Let g be a Lie algebra. When the dynamical form fields admit the inter- 
pretation of being components of a (higher) connection on a (higher) principal 
bundle defined by g over M, it becomes convenient to use g-valued forms as 
dynamical fields, that is, F — g <E) A k M. The Lie algebra naturally comes with 
the following bilinear operations: 

(commutator) [— ] : g ® g —> g, (91) 
(Killing form) (-) : g ® g -> R, (92) 

where [— ] is antisymmetric, while (— } is symmetric and it is non-degenerate for 
semi-simple g. They satisfy the following compatibility identity 

(a® [6<g>c]) = ([a (8) 6] 8 c), (93) 

where either side defines a trilinear, totally antisymmetric functional. We ex- 
tend slightly the notation for F-local horizontal forms and write £l h ' v (F, g) for 
the space of local g- valued variational (h, v)-forms; we also write f2 fe (M, g) = 
0(g>f2 fc (M). The operations d v , dh, *, [— ] and (— ) extend to local g-valued vari- 
ational forms simply by treating one or the other tensor factor trivially, while 
A is extended by acting as ® on the Lie algebra factor. Then it is easy to check 
that the expressions 

(Ai A *A 2 ) and A \p 2 A (jl 3 ]) (94) 
are totally symmetric in their arguments as long as /Ltj are forms of odd degrees. 
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